Lattice Glass Models 
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Motivated by the concept of geometrical frustration, we introduce a class of statistical mechanics 
lattice models for the glass transition. Monte Carlo simulations in three dimensions show that they 
display a dynamical glass transition which is very similar to that observed in other off-lattice systems 
and which does not depend on a specific dynamical rule. A mean-field study shows the existence of 
a discontinuous glass transition, in agreement with the numerical observations. 



Understanding the glass transition, and the glass 
phase, is one of the present major challenges in condensed 
matter physics. The experimental glass transition is re- 
lated to a dramatic dynamical slowing down in which the 
structural relaxation time changes of 14 orders of magni- 
tudes in a relatively small window of temperatures. As a 
matter of fact, the glass transition temperature T g is em- 
pirically defined as that where the structural relaxation 
time becomes of the order of an hour. In the last fifty 
years there have been many efforts to try to understand 
whether this phenomenon is just a dramatic crossover, it 
is related to a pure dynamical transition, or it is the sig- 
nature of a true thermodynamic glass transition, called 
the "ideal glass transition", which would take place at 
temperature Tk below T g but which is kinctically avoided 
|ij . The possibility of this last scenario in fragile glasses 
is supported by the closeness of the Vogel-Fulcher tem- 
perature (where the extrapolated relaxation time has a 
divergence) and of the Kauzmann temperature (where 
the extrapolated excess entropy of the supercooled liquid 
vanishes) M , and has been widely explored theoretically 
[Sir~0 ■ The analogy with mean field discontinuous spin 
glasses (6) has given recently a new boost to this line 
of research |tJ . These models can be solved analytically 
and display striking similarities with the experimental 
fragile glasses. However, they predict (because of their 
mean field character) a dynamic freezing transition at a 
temperature T c > T g , equal to the mode coupling tran- 
sition temperature, where real systems still have a finite 
relaxation time. This difference is generally attributed 
to the existence in finite dimensions of activated pro- 
cesses which would transform the 'mean field' dynamical 
transition into a crossover in such a way that the relax- 
ation time increase very fast in the temperature region 
Tk <T < T c and eventually diverges at Tk- 

In this paper we introduce new lattice models of 
glasses. We study them numerically in three dimensions, 
and show (1) that they display a dynamic glass transition 
similar to the one seen in the simulation of glass formers 
like binary Lennard- Jones ||, (2) that their mean field 
solution, obtained using the Bethe approximation, pre- 
dicts the same physical scenario found for mean field dis- 
continuous spin glasses (in particular a dynamic freezing 
transition preceding the "ideal glass transition"). 

Our lattice glass models are defined as follows. On 
each node of the lattice (e.g. a cubic lattice), there can 



be or 1 particle, but these occupations are restricted 
by a hard 'density-constraint': a particle cannot have 
more than I among its 6 neighboring sites occupied. One 
possible interpretation of this model is a coarse grained 
version of a usual (off-lattice) hard sphere system taking 
into account the effect of geometric frustration ||. One 
site of our lattice is then characterizing the density of 
spheres in a local cell of space, of volume of the order 
of an icosahedron built with the original spheres. The 
presence of a 'particle' in our model means that the lo- 
cal arrangement of spheres in this cell is very dense, as 
happens when icosahedral order sets in. The absence of 
the 'particle' corresponds to a less dense arrangement of 
spheres. Geometric frustration means that it is not pos- 
sible to fit together the high density icosahedral struc- 
tures to fill the space: this is taken into account by the 
'density-constraint ' . 

Notice the important difference with the Kob- 
Andersen model B Jl^] . In that case the jamming is forced 
by a dynamical rule: a particle which violates the density 
constraint is blocked. In our case, the model is defined 
thermodynamically: configurations violating the density 
constraint are forbidden. The thermodynamics defini- 
tion has two advantages which were absent in previous 
lattice models : 1) The existence of a dynamical 

phase transition, and the value of the density at which it 
takes place, do not depend on the type of local dynam- 
ics which is used, e.g. whether particles just hop on the 
lattice, or whether they are exchanged grand-canonically 
with a reservoir. 2) One can perform some analytic stud- 
ies of the thermodynamics, relate them to the dynamical 
observations, and address key issues like the possible ex- 
istence of an ideal glass transition. 

We have run some numerical simulations of the lattice 
glass models using two algorithms. The first one (CA) 
is a simple Monte-Carlo simulation at fixed density in 
the canonical ensemble, where a randomly chosen parti- 
cle can hop to a neighboring site if the density-constraint 
is satisfied. In order to find an acceptable initial con- 
figuration, we prepare the system through an annealing 
procedure, in which the constraint is soft: a particle with 
r neighbors has an energy E = (r — 1)6 {r — I) (8 is Heav- 
iside's step function), and we simulate the system with 
Metropolis algorithm at decreasing values of the temper- 
ature T. When a zero energy configuration is found, we 
turn to the hard density-constraint case (i.e. T = 0) and 
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start our canonical run. At equilibrium, the chemical 
potential is measured from p = hi(p/p), where p is the 
fraction of sites in which it is possible to add a new par- 
ticle and p is the density. The second type of simulation 
(GCA) uses the grand canonical ensemble. We have a 
reservoir with chemical potential p which is coupled to 
each lattice site and can create or destroy particles. 
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FIG. 1. The inverse of the density is plotted versus the 
chemical potential of the A particles for the mixture m\z on 
a 3-d cubic lattice of size 15 3 . Circles are obtained from the 
canonical simulation CA with 10 3 Monte Carlo Steps per Par- 
ticle (MCSP). Triangles, squares and diamonds are obtained 
from the GCA with, respectively, 10 3 , 3.10 3 , 10 5 MCSP at a 
fixed increasing rate of the chemical potential between p — 
and p = 8. In the last case, the transition to the crystalline 
state is clearly visible. Inset: The diffusion coefficient for 
the two types of particles, in the CA simulation of the same 
mixture, is plotted in log-log scale versus p c — p, showing the 
dynamical transition at p = p c . 

As expected, the GCA simulation reaches the equilib- 
rium faster. The obtained results look qualitatively very 
similar to the ones on Lennard- Jones systems For 
a one component fluid, we find for all £ > 1 that upon 
increasing p the system has a first order phase transi- 
tion towards a crystal, identified by a discontinuity of the 
density, and the presence of Bragg peaks in the diffrac- 
tion pattern. In order to study the glass transition, we 
have considered lattice glass binary mixtures, for which 
the tendency to crystallization is reduced. In these mix- 
tures, one obtains a reproducible "supercooled liquid" 
which exhibits, when the density increases, a dynami- 
cal glass transition. We present the results obtained for 
a mixture containing 30% of particles A with density- 
constraint I a = 1, and 70% of particles B with density- 
constraint Ib = 3 (called mi3 mixture in the following). 
Fig.|l| shows the density as a function of the chemical po- 
tential, measured in the two algorithms. One sees a clear 
saturation which takes place at a density p c ~ 0.565: 
none of these two local algorithms can reach a density 
higher than p c . The Inset of Fig shows the diffusion 
coefficient Da (resp. Db) of the particles measured in 



the CA simulation. Their decrease with increasing p is 
well fitted by a function vanishing at p c with a power 
law: 

Da,b{p) - CaAPc - p) aA ' B 0( P c - p) (1) 

with exponents of order oia ~ 4.2 and as ~ 2.3. 

This vanishing is generally taken as the numerical sig- 
nature of the dynamical phase transition p| , ^2[ . When 
the chemical potential is quenched above the critical 
value (corresponding to the dynamical transition) the 
systems remains out of equilibrium and has an aging be- 
havior which we have seen measuring the two time corre- 
lations of occupation numbers of a site. On the contrary, 
if we increase the chemical potential very slowly, we can 
observe the transition to a crystalline phase even in this 
case of mixtures as shown in Fig. [I]. It is quite possible 
that the 77743 mixture will eventually phase separate on 
very large times, but we have not seen any such effect on 
the time scales of our simulations. Whether there exist 
mixtures where the true equilibrium state at large p is a 
glass, with a thermodynamic phase transition from the 
liquid to this glass, is under current study fl3[| . Let us 
finally remark that the existence of a metastable glassy 
phase does not depend on the type of local dynamics we 
use (GCA or CA), only the time to nucleate the crystal 
does (contrary to the Kob- Andersen model). 

We now present a mean field theory of the finite di- 
mensional lattice glass models. We focus on the Bethe 
lattice version of these models |ll| , whose underlying lat- 
tice structure is a random graph with fixed connectivity: 
every vertex has exactly k + 1 neighbors, but the graph 
is otherwise random (to study the three dimensional case 
one takes k = 5). Locally (on finite lenghtscales), such 
a graph has the structure of a Cayley tree with a fixed 
branching ratio, but it also has loops of typical size In TV. 
The presence of these loops is crucial to induce the ge- 
ometric frustration, but the local tree structure allows 
for an analytic solution of the model. For the sake of 
clarity, we shall describe this solution in the simple case 
where there is only one type of particle, with a density- 
constraint given by the integer £• we have extended these 
computations to mixtures fllq] and we shall present the 
results for the 77713 mixture. 

Since the underlying lattice structure is locally tree- 
like, one can write iterative equations on the local proba- 
bility measure. More precisely, let us analyze one branch 
of the tree ending on site i. We denote by j 6 {1, k} 

(i) 

all the neighbors of i. We call Zq the partition function 
of this branch, restricted to configurations where site i 
is empty, and zfl the partition function restricted to 
configurations where site i is occupied, and has p neigh- 
bors occupied. Defining Sq = J2l=o %i,p> we nnc ^ the 
recursion relations: 

4 i) =n(^ ) +s«) 



2 



*2 = eM E 
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It is convenient to introduce on any site m the local fields 



/i m = ln(s< m V^ m) ) and a m = In 
in terms of which the iteration reads: 



n 1 + e h i 



at 



(3) 



where 



a 9 = E 

p=0 



l<jl<...<jp<fc 



(4) 



and the sum over ji < ... < j p is defined to take the value 
1 when p is zero. 

From these fields, one can obtain jl6| the grand canoni- 
cal potential A — — log(Z) as a sum of sites and links con- 

tribution A = -^A^ +E< M -> A% The 

contribution from the bond < i, i' > is given by the local 
fields obtained from the two branches arriving on i (in 
absence of i') and on i' (in absence of i) as: 



1 + e hi + e h ' 



ex P (-4?) 
The contribution from site i reads 



exp 



(-AS") 



= 1 



hi-\-ai~\-h i / + a^. 



(5) 



(6) 



where is the total field on site i, given by formula (g) 
where k is changed to k + 1. The shift in grand potential 
when merging the k branches ending at points j = X, k 
onto the node i is given by 



exp (-A'A = l + e h 



(7) 




FIG. 2. The entropy per site (continuous line) and the total 
density (dashed line) predicted by the liquid solution are plot- 
ted as a function of the chemical potential of the A particles 
for the mixture mi3 on a k = 5 Bethe lattice. 



Starting from the previous eqs., we find at low density 
a liquid phase characterized by a homogeneous solution 
hi = h,cii = a, where h, a verify two simple self-consistent 
equations. Given this solution one can easily obtain all 
the thermodynamic quantities, and in particular the den- 
sity p and the entropy per lattice site S = A — pji. For 
every k > 2 and t < k, (and for generic mixtures), the 
entropy becomes negative when the chemical potential (i 
becomes larger than a certain critical value p s =o- There- 
fore a thermodynamic phase transition takes place at a 
chemical potential /i < /J, s —o, as shown in Fig. |2|. 

To gain some further insight on the thermodynamic 
transition we study the stability of the liquid phase, i.e. 
we analyze all the generalized susceptibilities: 
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where U{ € {0, 1} is the occupation number of site i and < 
niUj > c = Gd denotes the connected correlation function 
between points i and j at a distance d. Since there is 
only one finite path (with probability one) connecting 
two points at a finite distance d , the computation of Gd 
(and therefore the stability analysis) can be reduced to a 
one dimensional problem which can be solved by transfer 
matrix technique [^p]| . The divergence of Si signals an 
instability toward a cyclic solution of the iteration eqs. 
fUSl , which corresponds to the crystal as found for I = 
by Runnels |1S| ]. Here, we focus on the next instability 
related to the divergence of the "glass susceptibility" 2 2 . 
Depending on the value oik, I (and the type of mixture), 
we can encounter two types of situations. The first one is 
when the glass susceptibility remains always finite (this is 
found for k = 1, 2, 3 and every I > 0), the glass transition 
is then discontinuous. Instead, when the susceptibility 
S2 diverges at /1 = /12, either there is a continuous glass 
phase transition at /i = U2, (this is found for £ = 0, 
which is nothing but the vertex covering model studied 
in |]l9| ) , or a discontinuous phase transition (without any 
divergence of S2) takes place at u c < /i2- This happens 
for sure whenever ^ s= o < M2 (this is found e.g. in the 
k = 5,1 = 3 case). The discontinuous transition case, 
on which we focus in the following, is the one relevant 
for glasses, contrary to the continuous transition which 
is typical of spin glasses 0. 

The high-density glassy phase can be analyzed in the 
Bethe approximation, taking into account the existence 
of many different local minima of the TAP free energy 
(called pure states in the following, see (^0|). In this case 
the fields hi and fluctuate not only from site to site 
but also from pure state to pure state. One defines for 
each site a probability distribution Ri(h, a) that the fields 
hi, a,i equal h, a for a randomly chosen pure state. In our 
case we have verified that Ri does not fluctuate from site 
to site, and the analysis of the high-density glassy phase 
reduces to obtaining a single function R(h, a). Using the 
cavity or the replica method [jl6| we find that this func- 
tion satisfies the self-consistent equation: 
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(l^-fe^)-" 1 = / 6 [dhjdajRihj^j)] S (h - hi) 6 (a - a t ) 

(9) 

where Af is a normalization constant, /ij, Oj are the local 
fields on site i obtained when merging k branches which 
carry the fields {hj, cij}, see eq. (||), and m is a Lagrange 
multiplier which fixes the value of the free energy den- 
sity of the pure states giving rise to R(h, a) Q . We have 
solved this eq. numerically using the algorithm of |jl5| . 
For some choices of k and t, we find a scenario identical 
to the one of discontinuous spin glasses, with two transi- 
tions: (1) a dynamical transition at pd, pd, (2) an equilib- 
rium glass transition, due to an entropy crisis a la Kauz- 
mann, at a certain density p eq > pd and chemical poten- 
tial p eq > fid- Precisely, when we increase p starting from 
the liquid phase (where Ru q (h, a) = §(hu q —h)8(aii q —a)), 
we first encounter at ju = /Xj a non trivial solution 
R^iassih, a) of (|9|) which appears discontinuously, signal- 
ing the existence of many pure states. The static (equi- 
librium) transition appears at a higher chemical potential 
when these new solutions dominate the thermodynamics 
JHi].^. For the m.13 mixture on a k — 5 lattice, we find 
Pd — 0.58, which is surprisingly close to the 3D value. 
The discontinuous character of the transition is assured 
by the fact that the difference q\ — qq is finite and positive 
at the transition, where the overlaps go and qi are defined 
as usual |2(|: q = (1/N) £\ < m > a < n* >p /N and 
qi = (l/N)J2i < n i >q /N and the indices a, {3 de- 
note two pure states randomly chosen according to their 
Boltzmann weights. Moreover we have checked, by run- 
ning a CA simulation on the Bethe lattice [|l6|, that the 
dynamical transition takes place at the density where 
the diffusion coefficients vanish. A detailed study of the 
thermodynamics of the lattice glass models, and of their 
equilibrium transition, will be presented elsewhere Jl| . 

In this letter we have introduced a new class of lat- 
tice models of glasses. Their 3D numerical simulations 
display a phenomenology very similar to the one of the 
Lennard-Jones systems, whereas their solution in the 
Bethe approximation exhibits a discontinuous glass tran- 
sition. These models allow to bridge the gap between 
the phenomenology of fragile glasses and their disordered 
spin glass analogs. They should allow to study key is- 
sues like the existence or not of a true thermodynamic 
phase transition, and the finite dimensional counterpart 
of the mean field dynamical transition, because (1) they 
are among the simplest systems to exhibit a glass transi- 
tion in 3D, (2) this transition is not linked to a specific 
dynamical rule, but it is present in any local dynamics, 
(3) contrary to off-lattice models, they can be solved in 
the Bethe approximation (i.e. in the limit of infinite di- 
mension) , (4) within such an approximation scheme they 
present a discontinuous spin glass transition. 
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